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H ■ 1 Abstract 
. ^ 

2 In this paper wc study a facility location problem in the plane in which a single point (facility) 

3 and a rapid transit line (highway) are simultaneously located in order to minimize the total travel 

4 time of the clients to the facility, using the Li or Manhattan metric. The rapid transit line is 

5 represented by a line segment with fixed length and arbitrary orientation. The highway is an 

6 alternative transportation system that can be used by the clients to reduce their travel time to 
, 7 the facility. This problem was introduced by Espejo and Rodn'guez-Chfa in [8]. They gave both 
' 8 a characterization of the optimal solutions and an algorithm running in 0{n^ logn) time, where 

^9 n represents the number of clients. In this paper we show that Espejo and Rodrfguez-Chfa's 

O ' 10 algorithm does not always work correctly. At the same time, we provide a proper characterization 

11 of the solutions and give an algorithm solving the problem in 0{n^) time. 

m 

• 12 Keywords: Geometric optimization; Facility location; Transportation; Time distance. 

!>■ . 

O ■ 13 1 Introduction 
^' 

' 14 Suppose that we have a set of clients represented as a set of points in the plane, and a service facility 

15 represented as a point to which all clients have to move. Every client can reach the facility directly 

16 or by using an alternative rapid transit line or highway, represented by a straight line segment of 

17 fixed length and arbitrary orientation, in order to reduce the travel time. Whenever a client moves 

18 directly to the facility, it moves at unit speed and the distance traveled is the Manhattan or Li 
cd • 19 distance to the facility. In the case where a client uses the highway, it travels the Li distance at 

20 unit speed to one endpoint of the highway, traverses the entire highway with a speed greater than 

21 one, and finally travels the Li distance from the other endpoint to the facility at unit speed. All 

22 clients traverse the highway at the same speed. Given the set of points representing the clients, 

23 the facility location problem consists in determining at the same time the facility point and the 

24 highway in order to minimize the total weighted travel time from the clients to the facility. The 
26 weighted travel time of a client is its travel time multiplied by a weight representing the intensity 

26 of its demand. This problem was introduced by Espejo and Rodriguez-Chia [8]. We refer to [8] and 

27 references therein to review both the state of the art and applications of this problem. 
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28 Geometric problems related to transportation networks have been recently considered in computa- 

29 tional geometry. Abellanas et. al. introduced the time metric model in [1]: Given an underlying 

30 metric, the user can travel at speed v when moving along a highway h or unit speed elsewhere. The 

31 particular case in which the underlying metric is the Li metric and all highways are axis-parallel 

32 segments of the same speed, is called the city metric [3]. The optimal positioning of transportation 

33 systems that minimize the maximum travel time among a set of points has been investigated in de- 

34 tail in recent papers [2l[6lll]. Other more general models are studied in |9]. The variant introduced 

35 by Espejo and Rodriguez-Chia aims to minimize the sum of the travel times (transportation cost) 

36 from the demand points to the new facility service, which has to be located simultaneously with a 

37 highway. The highway is used by a demand point whenever it saves time to reach the facility. 

Notation to formulate the problem is as follows. Let S be the set of n client points; / the service 
facility point; h the highway; I the length of h\ t and t' the endpoints of /i; and v >! the speed in 
which the points move along h. Let > be the weight (or demand) of a client point p. Given 
a point u of the plane, let x{u) and y{u) denote the x- and y— coordinates of u respectively. The 
distance or travel time (see Figure [T]) , between a point p and the service facility / is given by the 
function 

dt,t'{p, f) := min ( \\p - ||p - + - + - /||i, \\p-t'\\i + - + \\t - 

p • 



ft' 

/►--- 

Figure 1: The distance between a point p and the facility / using the highway. 

38 Then the problem can be formulated as follows: 

39 The Facility and Highvi^ay Location problem (FHL-problem): Given a set 5 of n 

40 points, a weight Wp> Q associated with each point p of S*, a fixed highway length ^ > 0, 

41 and a fixed speed > 1, locate a point (facility) / and a line segment (highway) h of 

42 length i with endpoints t and t' such that the function X^pg^ Wp ■ dt^t' (p, f) is minimized. 

43 Espejo and Rodriguez-Chia [8] studied the FHL-problem and gave the following characterization of 

44 the solutions. Consider the grid G defined by the set of all axis-parallel lines passing through the 

45 elements of S. They stated that there always exists an optimal highway having one endpoint at 

46 a vertex of G. Based on this, they proposed an 0(n'^ log n)-time algorithm to solve the problem. 

47 In this paper we show that the characterization given by Espejo and Rodriguez-Chia is not true in 

48 general, hence their algorithm does not always give the optimal solution. 

49 Addendum An anonymous referee pointed out that the authors of [8] published a corrigendum 

50 to their paper the 19th of January 2012, and that our result was not novel. In here we provide 

51 a chronological order of the events so that the reader can reach his/her own conclusions. The 
62 first version of this paper appeared the 5th of April 2011 on arXiv (and a preliminary version also 
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53 appeared in the proceedings of the Spanish Meeting on Computational Geometry the 27th of June 

54 2011). We contacted the authors of [8], and provided them with a copy of our paper, including the 

55 counterexample. Naturally, they were interested in our research, and wanted to know where had 

56 they done a mistake. The 29th of October 2011, the authors of [8] contacted us claiming that they 

57 had found the error in their paper. They provided us a write-up containing the corrected version 

58 of their proof, and suggested we combine our results. Given the difference in notation and the fact 

59 that this paper subsumes their result, we declined. From the conversation we can only deduce that 

60 the authors of [8] submitted their corrigendum sometime in early November 2011. 

61 As of now (16th March 2012), our paper is currently under supervision for journal publication, 

62 whereas the corrigendum has already appeared at COR. Although we would love if the submission, 

63 correction and publication process takes less than three months (as it appears to have happened 

64 with corrigendum at Computers and Operations -Researc/i journal), we understand that this is not 

65 possible in high-end journals. Regardless of our personal opinion of the actions of Espejo and 

66 Rodriguez-Chia, we believe that the date in which the result was found (and not published in a 

67 journal) is the relevant one. Thus, we claim that our paper is the first one to claim the error of [8]. 

68 On a side note, we note that the corrigendum of Espejo and Rodriguez-Chia is also wrong, since 

69 they claim that our characterization is weaker. They specifically say that "The description given by 

70 [this paper] means an infinite many number of candidates to be one of the endpoints of an optimal 

71 segment". Although Lemma [2. II does not explicitly say so, the algorithm of Section [3] only considers 

72 O(n^) cases (in particular a finite amount). 

73 Paper Organization In Section [2] we first provide a proper characterization of the solutions. 

74 After that we give a counterexample to the Espejo and Rodrfguez-Chia's characterization. We 

75 provide a set of five points, all having weight equal to one, and prove that no optimal highway has 

76 one endpoint in a vertex of G. In Section [3] we present an improved algorithm running in O(n^) 

77 time that correctly solves the FHL-problem. Finally, in Section [5j we state our conclusions and 

78 proposal for further research. 



79 2 Properties of an optimal solution 

80 A primary observation (also stated in j8j) is that the service facility can be located at one of the 

81 endpoints of the rapid transit line. From now on, we assume f = t' throughout the paper. This 

82 assumption simplifies the distance from a point p € S to the facility to the following expression, 

dt{p, f) = min | \\p - ||p - ^| . 

83 Using this observation, the expression of our objective function to minimize is <!>(/, t) = X^pg^u^p • 

84 dt{p, /). We call this value the total transportation cost associated with / and t (or simply the cost 

85 of / and t). 

86 We say that a point p uses the highway if||p — + | < — /||i, and that p does not use it (or 

87 goes directly to the facility) otherwise. Given / and t, we call travel bisector of / and t (or bisector 

88 for short) as the set of points z such that ||z — /||i = — + see Figure [2j A geometrical 

89 description of such a bisector can be found in [8] , as the boundary of the so-called captation region. 

90 Lemma 2.1 There exists an optimal solution to the FHL-problem satisfying one of the next con- 
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Figure 2: The bisector of / and t. 

91 ditions: 

92 (a) One of the endpoints of the highway is a vertex of G. 

93 (6) One endpoint of the highway is on a horizontal line of G, and the other endpoint is on a 

94 vertical line of G. 

95 Proof. Let / and t be the endpoints of an optimal highway h and assume neither of conditions (a) 

96 and (6) is satisfied. Using local perturbation we will transform this solution into one that satisfies 

97 one of these conditions. Assume neither / nor t is on any vertical line of G. Let 5i > Q (resp. 

98 62 > 0) be the smallest value such that if we translate h with vector (— 5i,0) (resp. ((52,0)) then 

99 either one endpoint of h touches a vertical line of G or a demand point hits the bisector of / and 

100 t. Given e G [— (5i,(52], let fe, t^, and he be /, t, and h translated with vector (e, 0), respectively. 

101 It is easy to see that \dt^{p, fe) — dt{p, f)\ = \e\ for all points p. Given a real number x, let sgn(a;) 

102 denote the sign of x. We partition S into three sets Si, S2 and 53 as follows: 

51 = {peS\ sgn{dtAp, fe) - dt{p, /)) = sgn(e), Ve G [-61,62] \ {0}} 

52 = {p G 5 I sgn(di,(p, /,) - dt{p, /)) = - sgn(e), Ve G [-61,62] \ {0}} 

53 = {p e S [ sgn{dtAp, fe) - dt{p, f)) = -1, VeG [-5i,<52]\{0}} 

103 Observe that points of 53 are in the bisector of / and t; Si contains the demand points that travel 

104 rightwards to reach / directly or by using the highway, and S2 contains the points that travel 

105 leftwards. 

106 Theoretically, one could consider the case in which a point belongs to set S^i = {p £ S [ sgn(dt^ {p, fe) — 

107 dt{p, f)) = 1, Ve G [— (5i,(52] \ {0}}. Geometrically speaking, the points of this set are those that, 

108 when translating the highway in either directions, the distance between them and the entry point 

109 of the highway increases. This situation can only happen when the point is aligned with the entry 
no point. That is, point p G 6*4 if and only if either [i) p uses the highway to reach the facility and it is 

111 vertically aligned with t, or [ii] p walks to the facility and it is vertically aligned with /. However, 

112 by definition of 61 and 62, no point of S can belong to (or enter) S4, during the whole translation. 

113 By the linearity of the Li metric, whenever we translate the highway e units to the right (for some 

114 arbitrarily small e, < e < 61), the highway will be e units closer for points in 52 U S'i, but e units 

115 further away for points of Si . Analogously, the distance to the facility decreases for points in Si U ^3 

116 and increases for points of 52 when translating h leftwards. Let A'' = X^pGSi "^p ~ Spe52 ^^-"^ 

117 k = Wp. Thus, for any vector (e,0), e G [— (5i,(52] \ {0}, the change of the objective function 
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118 when we translate the highway with vector (e, 0) is equal to the following expression: 

<^{fe,te)-Hf,t) = ^Wp-dt^{pJs)-^Wp-dt{p,f) 

pes pes 

= e Wp — e — |e| Wp 

peSi peS2 peSi 

= Ne-k\e\ 

119 Since we initially assumed that the location of h is optimal, we must have both A'^ = A; = 

120 (otherwise translating h rightwards or leftwards would result in a decrease of the objective fuction). 

121 In particular, we can translate h in either direction so that the cost of the objective function is 

122 unchanged. 

123 More importantly, observe that the value of k must remain on the whole translation: if at some 

124 point it becomes positive we can find a translation from that point that reduces the cost of the 

125 objective function. In particular, the set 5*3 must remain empty during the whole translation. Any 

126 point that changes from set Si to 5*2 (or vice versa) must first enter S3. Since the latter set remains 

127 empty during the whole translation, no point can change between sets Si,S2, or S3 until either / 

128 or t is vertically aligned with a point of S. 

129 We perform the same operations on the y coordinates and obtain that one of the two endpoints is 

130 on a horizontal line of G, hence satisfying one of the two conditions of the Lemma. □ 

131 When the highway's length is equal to zero, the FHL-problem is the weighted 1-median prob- 

132 lem in metric Li [7], and in this case the item (a) of Lemma 12.11 holds. Espejo and Rodriguez- 

133 Chia ^ claimed that there always exists an optimal solution of the FHL-problem that satisfies 

134 Lemma [2. II (a). Unfortunately, this claim is not true in general and their algorithm may miss some 

135 highway locations; indeed, it may miss the optimal location and thus fail. We provide here one 

136 counterexample and the following result. 

137 Lemma 2.2 There exists a set of unweighted points in which no optimal solution to the FHL- 

138 problem satisfies Lemma \2.1\ (a). 

139 Proof. Consider the problem instance with five points whose coordinates are (—4,0), (—3,-1), 

140 (12, 8), (13, 5), and (13, 7), respectively (see Figured]). In the problem instance, we give unit weight 

141 to all points, and set the length h of the highway as ^ = VTSO ~ 13, 5. For simplicity in the 

142 calculations, we also set v = i, but any other large number works as well. The cost associated to 

143 the highway of endpoints / = (12,6) and t = (0,0) is 10 + 2£/v = 12. We claim that this location 

144 is better than any other solution with an endpoint at a vertex of G. 

145 If one endpoint of /i is a vertex of G in the line x = —3, then the other endpoint is located to the 

146 left of the line x = 11 because —3 + £ < 11. In that case we can translate h rightwards with vector 

147 (^, 0) improving the objective function. The same holds if one endpoint of /i is a vertex in the line 

148 X = —4. Similarly, if one endpoint is a vertex in the line x = 13, then we can translate h leftwards 

149 with vector (—^,0) and the objective function decreases. 

150 Consider now locating one of the highway endpoints at coordinates (12, 0) or (12, —1). Observe that 

151 the walking time (i.e., the traveling time when the highway is not used) from the points (—4, 0) and 

152 (—3, —1) takes at least 15 units of time, which is more than the cost associated with our solution. 

153 The same happens to the sum of the traveling times of the three other points. Hence, if / is located 

154 at one of the two vertices, the five points must use highway (otherwise the travel time is higher than 
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Figure 3: A counterexample to the algorithm of Espejo and Rodn'guez-Chia. 



155 our solution). Analogously, if t is located at grid points (12, 0) or (12, —1), no point of S will use the 

156 highway. In either case, the corresponding solution is at least as high as the sum of distances from 

157 all points of S to the geometric median, which is higher than the cost associated with our solution. 

Consider now the cases in which one of the endpoints has coordinates (12, yo) for some yo £ {5, 7, 8}. 
We start by showing that, in any of the three cases, the optimal position of the other endpoint of 
the highway (denoted by e) must lie on the line y = 0. Since the highway's length is equal to i, 
the possible positions of e lie both in circle a of radius i centered at (12, yo) and to the left of line 
X = 12. Observe that the clients that walk to e are points a = (—4,0) and b = (—3,-1), located 
always to the left of e. Hence, we are interested in minimizing the expression ||a — e||i + ||6 — e||i. Let 
a' ,b' € a denote respectively the closest points to a and b with the Li metric, which verify y{a') = 
and y{b') = —1. Observe that if y{e) > then ||a — a'||i < ||a — e||i and ||6 — a'||i < — e||i implying 

||a — a'||i + ||6 — a'lli < ||a — e||i + ||6 — e||i 

(see Figure m a)). Similarly, if y{e) < —1, then 

||a — + — 6'||i < ||a — e||i + ||6 — e||i. 

158 Therefore, e must satisfy —1 < y(e) < (see Figure Hlb)). In this case we have 

||a - e||i + ||6 - e||i = x{e) - x{a) + y{a) - y{e) + x{e) - x{b) + y{e) - y{b) 
= 2x{e) + 8 

159 Then ||o — e||i + ||6 — e||i is minimized when x{e) is minimum, and it happens when y{e) = 0. 

160 If yo = 8, then h can be translated downwards with vector (0, — |) and the value of the objective 

161 function decreases. Thus point (12, 8) is discarded. It remains to show that there is a solution 
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b ^ 
%/ / a) 

Figure 4: a = (—4, 0) and h = (—3, —1). When one endpoint of the highway has coordinates (12, 8), (12, 7), 
or (12, 5), the optimal position of the other endpoint e is on the line y = 0- 

162 better than the one having an endpoint at either (12,7) or (12,5), and the other endpoint on the 

163 line y = 0. Observe that if / and t belong to the lines y = and x = 12, respectively, then by 

164 exchanging / and t the value of the objective function reduces in £/v. Then consider the case where 

165 y{t) = and x(/) = 12. 

166 Let t = (0,0) and / = (12,6). Given a value e, let be the point with coordinates (e, 0) and fs 

167 be the point in the line x = 12 such that ^(/e) > and the Euclidean distance between and 

168 te is equal to i (see Figure [5]). Let [— (^i,(52], Si, 62 > 0, be the maximal-length interval such that 

169 5 < y{fe) < 7 for ah e G [-61,62]. Note 61 = \/T55 - 12 < 1 and <52 = 12 - \/l3T < 1. Then |e| < 1. 
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Figure 5: Definitions of fe and te- 

no 

171 The variation of the objective function's value when / and t are moved to fs and t^, respectively, 

172 is equal to 

gie) := <&(/„ t,) - <^>(/, t) 

= 2{x{te) -x{t))-{y{fe )-y{f)) 
= 2e- (^\/36 + 24e-e2 -e) . 

173 In the following we will show that V36 + 24e — < 6 + 2e, for all e S [—61, 62] \ {0}. In particular, 

174 we will have g{£) > (except when e = 0), implying that our highway location is optimal. First 

175 observe that Ae"^ + 24e + 36 = {2e + 6)^ > 36 + 246 - e^ . Since \e\ < 1 then 2e + 6 > and 

176 36 + 24e — > 0, which implies 2e + 6 > ^36 + 24e — e^. Thus 5(e) > and the highway with 
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177 endpoints / and t gives a better solution than that having an endpoint at (12,7) or (12,5). This 

178 completes the proof. □ 



179 In the next section we provide a correct algorithm that solves the problem in 0{n^) time. We 

180 assume general position, that is, there are no two points on a same line having slope in the set 

181 {— 1,0,1, oo}. 



182 3 The algorithm 

183 Lemma |2. II can be used to find an optimal solution to the FHL-problem. Although the method is 

184 quite similar for both cases in Lemma l2.H we address the two cases independently for the sake of 

185 clarity. By Vertex-FHL-problem we will denote the FHL-problem for the cases in which Lemma [2. II 

136 a) holds, and by Edge-FHL-problem the FHL-problem for the cases in which Lemma l2. li b) holds. 

137 In the next subsections we give an 0(n^ )-time algorithm for each variant of the problem. In both 

138 of them we assume w.l.o.g. that highway's length i is equal to one. 

139 In the following 9 will denote the positive angle of the highway with respect to the positive direction 

190 of the X-axis. For the sake of clarity, we will assume that 9 G [0, j]. When 9 belongs to the interval 

191 [kj, + 1)^], A; = 1, . . . , 7, both the Vertex- and Edge-FHL-problem can be solved in a similar way. 

192 Given a point u and an angle 9, let u{9) be the point with coordinates {x{u) + cos 9, y{u) + sinS). 

193 There exists an angle (j) € [0, j] such that the bisector of the endpoints / and t = f{9) has the shape 

194 in Figure [2] a) for all 9 G [0, (^), and has the shape in Figure [2] b) for all 9 G (0, j]. Such an angle (p 

195 verifies cos((/>) — sm{(j)) = K Furthermore, 4> = ^ arcsin(l — \) and (/) j unless v is infinite. Refer 

196 to [8j for a detailed description of this situation. 

197 Let Hj., Hj^, and Hx+y denote the point set S sorted according to the x-, y-, and (x + y)-order, 
193 respectively. 



199 3.1 Solving the Vertex-FHL-problem 

200 For each vertex u of G we can solve the problem subject to f = u or t = u. We show how to obtain 

201 a solution if / = u. The case where t = u can be solved analogously. 

202 Suppose w.l.o.g. that the vertex / = n is the origin of the coordinate system and the highway angle 

203 is 9, for 9 G [0, j]. Then t = u{9) = {cos9,sm9). and the distance dt{p, f) between a point p £ S 

204 and the facility u has the expression ci + C2Cos9 + cssin^, where ci,C2,C3 are constants satisfying 

205 C2, C3 G { — 1, 0, 1}. When 9 goes from to ^, this expression changes at the values of 9 such that: 



The point p switches from using the highway to going directly to the facility (or vice versa). 
We call these changes bisector events. A bisector event occurs when the bisector between the 
highway's endpoints u and u{9), contains p. At most two bisector events are obtained for each 
point p. 

The highway endpoint u{9) crosses the vertical or horizontal line passing through p. We call 
this event grid event. Again, each point of S generates at most two grid events. 

9 = <j). We call it the (p-event. 
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213 We refer the interested reader to [8] for a detailed description of the above event^. The cost of 

214 their algorithm is dominated by the time spent sorting the order in which events take place. In 

215 order to avoid this time, we use the following result: 

216 Lemma 3.1 After an 0{n\ogn)-time preprocessing, the angular order of all the events associated 

217 with a given vertex of G can he obtained in linear time. 

218 Proof. The preprocessing consists in computing II^;, 11^,, and Hx+y, which can be done in 0(n log n) 

219 time. Now, let u be a vertex of G. It is straightforward to see that they are 0{n) grid events and 

220 that we can obtain their angular order in linear time by using both 11^ and liy. Let us show how 

221 to obtain the bisector events in 0(n) time. 

222 The bisector of u and u{9) consists of two axis-aligned half-lines and a line segment with slope 

223 -1 connecting their endpoints (see Figure [2] and [8] for further details). Given a point p, when 9 

224 goes from to 7r/4 the bisector between u and u{9) passes through p at most twice, that is, when 

225 p belongs to one of the half-lines of the bisector and when p belongs to the line segment. If p 

226 belongs to the line segment of the bisector then the event is denoted by ap (see Figure [6] b)). If p 

227 belongs to the leftmost half-line of the bisector, which is always vertical, we denote that event by 

228 Pp (see Figure El a)). Otherwise, if p belongs to the rightmost half- line which can be either vertical 

229 or horizontal we denote that event by 7p (see Figure [6] c) and d) ) . Observe that if the rightmost 

230 half-line is vertical then 7^ < 0, otherwise 7j, > (p. Refer to [8] for a characterization to identify 

231 whether a point p & S generates a bisector event for some angle 6. 




Figure 6: The bisector events of p when 9 £ [0, j]. a) p belongs to the leftmost half-line of the bisector of u 
and u{6). b) p belongs to the segment. c,d) p belongs to the rightmost half-line of the bisector. 

232 Let Hi be the subsequence of Ilx+y containing all elements p such that Op € [0, j], II2 be the 

233 subsequence of 11^ containing all elements p such that /3p € [0,^], and II3 be the subsequence 

234 of Hx that contains all elements p such that y{p) < y{u) and 7p G [O,*/*], concatenated with the 

235 subsequence of Hy that contains all elements p such that x{p) > x{u) and jp € [4>, j]. Given a point 

236 p G S, the corresponding events of p in [0, j] can be found in constant time, thus Hi, 112, and II3 

237 can be built in linear time. 

^Although their events are very similar to the ones we described, the authors of [8] refer to them as projection and 
limit points. We prefer to use the term "event" , since "point" is reserved for the elements of S 
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238 The following statements are true for any point p £ S: 

239 (a) x{p) + y{p) = ^(cos ap + sin ap + ^) for all points p in Hi. 

240 (6) x{p) = ^(cos /3p — sin f3p + i) for all points p in 112. 

241 (c) x(p) = |(cos7p + sin7p + ^) for all points p in Ha such that •jp < (j). 

242 (d) = ^(— cos7p + sin7p + ^) for all points p in lis such that jp > (j). 

243 Let Fi (resp. T^) be the sequence obtained by replacing each element p in Hi (resp. 112, Hs) 

244 by Up (resp. /3p, 7p). Therefore, from statements (a) — (d) and the monotonicity of the functions 

245 cos 9 + sin 0, cos 9 — sin 0, and — cos 9 + sin 9 in the interval [0, j\ , we obtain that Fi, F2, and F3 are 

246 sorted sequences. Using a standard method for merging sorted lists, we can merge in linear time 

247 Fi, F2, F3, the grid events, and the i;^-event. Therefore, the angular order of all events associated 

248 with a vertex u can be obtained in 0(n) time and the result follows. □ 

249 Theorem 3.2 The Vertex-FHL-problem can be solved in 0{n'^) time. 

250 Proof. Let u be a vertex of G. Using Lemma 13.11 we obtain in linear time the angular order of 

251 the 0(n) events associated with u. The events induce a partition of [0, j] into maximal intervals. 

252 For each of those intervals, the objective function takes the form g{9) := ^{f,t) = ^{u,u{9)) = 

253 hi + 62 cos + 63 sin 0, where ^1,62)^3 are constants. This problem is of constant size in each 

254 subinterval and the minimum of g{9) can be found in 0(1) time. Furthermore, the expression of 

255 g{9) can be updated in constant time when 9 crosses an event point distinct of <j) when going from 

256 to ^. In the case where 9 crosses (p, g{9) can be updated in at most 0(n) time. Then the problem 

257 subject to f = u can be solved in linear time. The case in which t = u can be addressed in a similar 

258 way. It gives an overall O(n^) time complexity because G has O(n^) vertices. □ 

259 3.2 Solving the Edge-FHL-problem 

260 We now consider the case in which the optimal solution satisfies condition b) of Lemma [2.1[ Namely, 

261 we consider a horizontal line of G and each vertical line Cy of G. For every pair of such lines, 

262 we consider eight different sub-cases, depending on whether h is located above/below eh, right- 

263 wards/leftwards of e^, and f £ eh and t £ e^ (or vice versa). For a fixed sub-case, we parametrize 

264 the location of the highway by the angle 9 that the highway forms with Ch ■ As in the Vertex-FHL 

265 case, we assume that f £ Ch, t £ Cy, and 9 £ [0, j]. 

266 We implicitly redefine the coordinate system so that Ch and Cy intersect at the origin o. Let 9 £ [0, j] 

267 be the positive angle of the highway with respect to the positive direction of the a;-axis and f = xg, 

268 t = yg he the highway endpoints, see Figure [3 

269 First notice that, since we are again doing a continuous translation of h, the events that affect 

270 the value of the objective function are exactly the same as those that happen in the Vertex-FHL- 

271 problem: bisector-, grid- and (p- events. We start by showing that the equivalent of Lemma 13.11 also 

272 holds: 

273 Lemma 3.3 After an 0{n\ogn)-time preprocessing, the angular order of all the events associated 

274 with a pair of perpendicular lines of G can he obtained in linear time. 
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Figure 7: Solving tlie Edgc-FHL-problem. 



275 Proof. We can follow the arguments of Lemma l3.11 Firstly, we note that there are 0(n) grid events 

276 and their angular order can be obtained in linear time by using both JJ^ and Hy. 

277 Given a point p £ S, let the events Op, I3p, and 7p be defined as in the Vertex-FHL case. Refer to 

278 Figure [6l Let Hi be the subsequence of ILx+y containing all elements p such that ap G [0, 112 be 

279 the subsequence of Hx containing all elements p such that f3p € [0, |], and Hs be the subsequence 

280 of Hx that contains all elements p such that y{p) < y{o) and 7p S [0, concatenated with the 

281 subsequence of that contains all elements p such that x{p) > x{o) and 7p G [</>, j]. Note that 

282 Hi, II2, and can be built in linear time. 

283 Given a point p G S, the following statements are true: 

284 (a) x{p) + y{p) = \{— cos Up + sin ap + ^) for all points p in Hi. 

285 (6) x{p) = \{— cos /3p — sin /3p + i) for all points p in 112. 

286 (c) x{p) = ^(— cos7p + sin7p + i) for all points p in such that 7p < (f). 

287 [d) y{p) = \{— cos7p + sin7p + for all points p in such that 7p > (j). 

288 Let Pi (resp. r2, T^) be the sequence obtained by replacing each element p in Hi (resp. n2, 

289 Ha) by Op (resp. /3p, 7fe). Therefore, by using similar arguments to those used in Lemma l3.ll the 

290 angular order of all events can be obtained in 0{n) time, once the lists n^;, 11^ and ^x+y have been 

291 precomputed. □ 

292 Consider now a small interval [^1,^2] in which no event occurs. Observe that, after the coordinate 

293 system redefinition, we have f = xq = (—cos 0,0), and t = yg = (O,£sin0). Let p € S he a point 

294 that uses the highway to reach the facility; since only the y-coordinate of t changes, its distance to 

295 / can be expressed as ci it sinO for some ci > 0. Analogously, if p walks to /, its distance is of the 

296 form ci lb cos 6 for some ci > 0. That is, the distance between a point of S and / in any interval is 

297 of the form ci + C2 sin 6 + 03 cos 9 for some constants ci > and 02,03 G { — 1,0,1}. 

298 Theorem 3.4 The Edge-FHL-problem can be solved in 0{n^) time. 

299 Proof. We can use a method similar to the one used in the Vertex-FHL-problem. Let eh be a 

300 horizontal line of G and e„ be a vertical line of G. 

301 Using Lemma Is. 3| we obtain in linear time the angular order of the 0{n) events associated with eu 

302 and The events induce a partition of [0, ^] into maximal intervals. For each of those intervals the 

303 objective function has the form g{9) := $(/, t) = ^{xq, yo) = 61 + 62 cos + 63 sin 6*, where 61 > 0, 
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Figure [8] 




V 


Cost 


Ratio 


a) 





- 


49 


1 


b) 


1 


2 


46 


0.93 






4 


45.5 


0.92 






10^ 


45 


0.91 


c) 


7.07 


2 


34.07 


0.7 






4 


30.54 


0.62 






10^ 


27 


0.55 


d) 


13.41 


2 


27.41 


0.56 






4 


20.71 


0.42 






10^ 


14 


0.29 


e) 


16.55 


2 


9 


0.18 






4 


8.5 


0.17 






10^ 


8 


0.16 



Table 1: Total transportation cost as a function of the highway's length and speed. The last column shows 
how much does the highway improve the total transportation cost (compared to the case in which only a 
facility is located) 



304 and 62,63 € Z are constants. This problem has constant size, hence the minimum of g{9) can be 

305 found in 0(1) time. Furthermore, the expression of g{9) can be updated in constant time when 9 

306 crosses an event point distinct of (j) when it goes from to ^. In the case where 9 crosses (j), g{9) 

307 can be updated in at most 0{n) time. Then the problem subject to / G e/,, and t £ can be solved 

308 in linear time. It gives an overall O(n^) time complexity because G has O(n^) pairs consisting of a 

309 horizontal and a vertical line. □ 



310 4 Experimental results 

311 Similar to [8] , we explore examples of solutions to the FHL-problem for different values of the length 

312 of the line segment. The problem instance is given by the unweighted points with coordinates 

313 (—4,0), (—3,-1), (12,8), (13,5), and (13,7) as in Lemma [2^2] and we consider locating a highway 

314 for different values of length and speed. Given a fixed value of speed, say v = 2, Figure [8] shows the 

315 location of the optimal highways for some values of L Note that the case i = is the Fermat- Weber 

316 problem for the Li-metric. The highway's length and the associated total transportation cost for 

317 each of these solutions can be seen in Table HI The optimal solution for each of the cases (and its 

318 associated cost) has been obtained with the help of a computer. 

319 Observe that, for some values of i, the optimal solution satisfies condition (a) of Lemma [2. H but in 

320 other situations condition (6) is satisfied instead (see Figure [8] d), where the highway's length has 

321 been set to 13.41). Experimentally we observed that increasing the highway's length decreases the 

322 total transportation cost until £ = \/305, in which a total cost of 5 + 2i/v is obtained (see Figure [8] 

323 e)). Afterwards the cost gradually increases until we locate a highway so long that no point of S 

324 uses it to reach /. We also note that for this demand point set the highway's speed has a small 

325 impact on the optimal solution. Indeed, increasing the highway's speed changes the total cost but 

326 the location of the highway in the above instance is unaffected by the highway's speed (provided 

327 that V > 1). The fourth column in Table H] gives the small variation of the total cost with respect 

328 to the speed. This suggests the following open problem: given an instance of the FHL-problem, can 

329 we efficiently compute the highway's length that minimizes the total transportation cost? 



12 




Figure 8: Solution of the same instance of the FHL-problem for different values of (!.. The optimal highway 
is depicted in red (and the endpoint containing / as a cross). The exact highway's length and the associated 
total transportation cost can be seen in Tabic 01 



330 5 Concluding remarks 



331 As further research, it would be worth studying the same problem in other metrics or using different 

332 optimization criteria. Another interesting variant would be to consider the problem when the length 

333 of the highway is not given in advance and it is a variable in the problem. Additionally, we could 

334 consider a similar distance model in which the clients can enter and exit the highway at any point 

335 (called freeway in [5]). 

336 Motivated from the experimental results of Section (H we can deduce that the highway's length 

337 has a strong impact on the optimal solution. As one would expect, when the highway's length is 

338 small, the total cost barely changes. We obtain a similar effect when the highway to locate is very 

339 long, since traveling to the opposite endpoint takes more time than walking directly to the facility. 

340 Hence, it would be interesting to consider a variation of the problem in which we can also adjust 

341 the highway's length. Specially, one would like to find a balance between the cost of constructing a 

342 longer highway and the improvement in the total transportation cost. 
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